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Studies have shown that standard lectures and instructional laboratory experiments are not effec- 
tive at teaching interference and diffraction. In response, the author created an interactive computer 
program that simulates interference and diffraction effects using the Finite Difference Time Domain 
method. The software allows students to easily control, visualize, and quantitatively measure the 
effects. Students collected data from simulations as part of their laboratory exercise, and they 
performed well on a subsequent quiz — showing promise for this approach. 



I. INTRODUCTION 

Often, interference and diffraction is taught by showing 
the mathematics in lecture and then performing laser in- 
terference and diffraction experiments in an instructional 
laboratory. However, studies show this teaching method 
is lacking. 

Consider the following problem (see Figure [T]) , which 
should be straightforward for a student who understands 
the concepts. Two point sources, 2.5 wavelengths (A) 
apart, generate waves in phase. Is there constructive 
interference, destructive interference, or neither at points 
A, B, and C? 

A University of Washington study asked this and other 
questions to « 1200 undergraduates in their standard, 
calculus-based, introductory physics course. Only « 35% 
answered correctly for both points A and B, and only 
~ 5% answered correctly for point C. Graduate students 
also performed poorly; only « 25% answered correctly for 
point CP These mistakes were primarily due to funda- 
mental misunderstandings; post-quiz student interviews 
included responses like, "I suppose that 2. 5 A [the slit 
separation] is small compared to 400 A and 300 A, so the 
sources here act like a single source. '^Moreover, the typ- 
ical instruction method did not help; the scores for point 
C before and after lecture and laboratory were basically 
unchanged 

Part of the problem is that light's wave nature is only 
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Figure 1: The diagram for the interference quiz. 



visible indirectly — though effects like interference and 
diffraction. Light's peaks and valleys cannot be viewed 
directly. Hands-on experience seeing and controlling the 
waves' interactions is useful for student understanding, so 
a better wave medium is needed. Ripple tanks provide 
this, but there are significant trade-offs between easy of 
use, cost, and measurement accuracy. 

For example, Wosilait et al. provided students with 
simple ripple tanks — pans, dowels, and sponges. They 
are cheap and easy to use but do not allow quantitative 
measurements. Still, they proved effective as part of an 
intensive tutorial system that raised scores for point C to 

* 45%.n 

Measurements can be performed with a ripple tank — 
at the price of increased cost and complexity. They re- 
quire a strobe device and a consistent wave source — not 
a student rolling a dowel back and forth. The strobe 
can either be a lightP (easier but more expensive) or a 
spinning, slotted disk 4 (cheaper but more complicated). 
Either makes the waves appear frozen, making their mea- 
surement easier. 

Computer simulations can also display waves, and they 
have potential advantages over ripple tanks. Simulations 
can allow finer control and more accurate measurements. 
They are available for free. They can be used outside 
of laboratory. And they allow better visualization by 
avoiding unwanted reflections, being easily paused, and 
having superior contrast. 

Dozens of simulations are available (search |http : | 
//www.merlot . org/). Some display what would be 
seen in a laser experiment — without showing the waves 
interacting.^ Others are non-interactive animations.^ 
Some have interactive displays!^ A few allow measure- 
ment of the instantaneous wave height/field strengthP 

Those simulations are primarily used for qualitative 
understanding. However, I wanted a simulation that 
could answer questions like the aforementioned one — 
albeit on a smaller scale — and could thereby test inter- 
ference conditions. This requires easily measured am- 
plitudes at any user desired grid point. No simulations 
combined that with interactive, easy to control, animated 
waves. 

Similar simulations have been proposed, but their 
use in instructional laboratories has not been reported. 
Frances et al. proposed using a Finite Difference Time 
Domain (FDTD) electromagnetics simulation to show in- 
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terference and diffraction effects, but they have only re- 
ported its use in lecture demonstrationsP Werley et al. 
also suggested using FDTD simulations in laboratories, 
but they instead used pre-recorded videos of actual prop- 
agating radiationP3 

So, I wrote an FDTD simulation that allows easy and 
accurate measurements, and I constructed a laboratory 
exercise that uses the simulation for quantitative mea- 
surements, not just qualitative understanding. 



II. THE PROGRAM 

The FDTD technique solves differential equations by 
discretizing them in both space and time. It has proved 
popular and effective for simulating electromagnetics, 
and there are many fine references on the techniqueP"^ 
So, the equations are not reproduced here; 1 this section 
summarizes aspects of the simulation and interface rele- 
vant to its use. 

Figure [2] shows the program's interface, which has 
five plots. The two large plots show E z (upper plot) 
and E z rms (lower plot). For both, black is the small- 
est value and white is the largest value, with shades of 
gray in between. The three smaller plots show E z and 
±\^2E z rms — an envelope for E z — along the horizontal 
and vertical dashed lines through the two larger plots. 
Those lines can be moved with the keyboard or mouse, 
and x, y, E Z: and E z rms at those lines' intersection is 
displayed in the center right area between the two verti- 
cal plots. Knowing E z rms a t that point allows users to 
home in on extrema. 

A plane wave — with a wavelength of 20 grid cells — 
enters from the left. It is not simulated but calculated 
analytically. At the start of the simulation, the wave's 
magnitude is ramped up gradually to avoid potentially 
unstable high frequency components. 

The barrier — the red line visible in both large plots — 
is a perfect conductor, and the openings in the barrier 
are hard sources that inject the incoming wave in to 
the FDTD domain — the area to the right of the barrier. 
Openings in the barrier can be added, removed, and mod- 
ified using the barrier control frame — at the right of the 
interface. 

Split-field perfectly matched layers^IH terminate the 
other three sides of the FDTD domain. These boundaries 
reduce reflections to imperceptible levels, effectively giv- 
ing the simulation open boundaries. 

When the simulation starts or the barrier is modified, 
a timer appears over the E z rms plot — counting down 
until steady state is reached. Afterwards, E z rms is reset 
to remove transients, and another countdown appears for 




Figure 2: The single narrow slit simulation. 



one wave time period. The steady state E z rms is calcu- 
lated by averaging over that time. 

Among its other features, the simulation also has a 
fast forward mode, which saves time by not updating the 
plots. In that mode, the simulation runs until the current 
countdown is done. 

The laboratory's computers — running Windows 7 with 
Intel Core 2 Duo processors — take ~ 55ms per timestep, 
resulting in a smooth animation. 

The software is written in Python using NumPy 
for the calculations, Tklnter for the interface, and 
the Python Imaging Library for the plots. Those 
libraries are available for Windows, OS X, and 
Linux. Executables, source code, and program infor- 
mation are available at |http : //lnmaurer . github . comTl 
Interference- Inference- Interface/; the program's 
source code available under the GNU Public License ver- 
sion 2. 



III. THE SIMULATIONS 

The laboratory was tailored for a particular class — 
aimed at future physics majors, but tweaking the exer- 
cises for other courses should be straightforward. Besides 
simulations, the laboratory also included pen and paper 
work and short laser /slit experiments. 2 However, four 
simulations are at the laboratory's heart: of narrow sin- 
gle, double, and triple slits and a wide single slit. 



1 For reference, the program simulates a TM Z wave using the stan- 
dard Yee lattice and second-order-accurate finite difference ap- 
proximations. 



The worksheet is available 

https: / /github. com/lnmaurer/Interference-Diffraction- 
Worksheet 



at 
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A. Narrow Single Slit 

See Figure [2] This simulation provides a baseline for 
comparison. It shows that a single narrow opening is 
not responsible for the interference patterns seen in the 
following simulations. 

The course's lectures included a mathematical descrip- 
tion of point sources in 3D, so this simulation introduced 
students to the asymmetric sources used in the labora- 
tory. To acquaint students with the simulation's controls, 
I had them roughly measure how fast the amplitude de- 
creased with distance from the opening; because this is 
a 2D simulation, the amplitude falls off less quickly than 
in 3D. 



B. Narrow Double Slits 

See Figure[3j This is the key simulation; it lets students 
discover the conditions for constructive and destructive 
interference. Students were asked to find four maxima 
and four minima of E z rms — m either x or y and in 
the right half of the domain, calculate \ d °~ dl \ for each — 
where do and d\ are the distances from the extrema to 
the openings, and then find the pattern in those numbers. 

In principle, the students already knew the conditions, 
but they seemed new to many. Furthermore, most stu- 
dents seemed to be comprehending the conditions for the 
first time. 3 Pen and paper work followed to reinforce 



3 Discretization in the program can slightly alter the conditions. 



their understanding. This included drawing a couple 
waves that re sulted in those conditions (using a simpli- 
fied version o f 15 ! 16 " ) and mathematically verifying that 
the conditions lead to waves in or out of phase. 



C. Narrow Triple Slits 

See Figure [4j The additional slit makes students 
closely consider the logic behind the extrema conditions. 

Here, the students were asked to find points that were 
simultaneously extrema in both x and y. This is straight- 
forward with the simulation, but — since this requires fine 
control in both x and y — this would be difficult using the 
standard experimental setup. 

The condition for constructive interference still holds. 
To get a maximum, all pairs should be maximized. 

However, the minimization condition is more 
complicated.^ The two slit condition caused the 
two waves to roughly cancel. That would leave the third 
wave undiminished. Additionally, it is mathematically 
impossible for all differences in the distances to be a half 
integer number of wavelengths. 

This is an instructive example that is missing from 
laboratories that do not use simulations for quantitative 
measurements. 



However, \do — d±\ was always within one. For example, the 
point where the dashed lines crosses is a maximum in Figure [3] 
and — x dl ■ = l 235 ^ 296 l — |l _ 3 Q5 These small errors were 
not problematic for students. 
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Figure 5: The wide single slit simulation with the dashed 
lines intersecting at the first minimum. It occurs at 6 — 
arctan ( H§E^ ) — 0.213, whereas the far field limit has 
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Figure 6: Figure for the second quiz question. 



B — students in the old University of Washington system 
scored 35% P The quiz included another question from 
the University of Washington study. See Figure [6] What 
would be observed at the far away point — first with just 
the left two slits and then with all three slits? 80% got 
both parts right. 



D. Wide Single Slit 



CONCLUSION 



See Figure |5j The right side of the domain gets close to 
the far field limit. While the domain is too small to truly 
display far field effects — a limitation common to many 
ripple tanks, the simulation successfully shows diffraction 
effects arising from waves. The simulation could be used 
to investigate near field effects — something difficult with 
normal instructional laboratory equipment. 



IV. RESULTS 



Table I: Scores for point C. The first three columns of results 
are from Wosilait et al. and are rounded to the nearest 5%. 1 



Interference and diffraction are tricky subjects to 
teach, but this approach shows promise in making them 
more accessible to students. The simulation allows them 
to easily control, visualize, and quantitatively measure 
interference and diffraction effects. That makes the soft- 
ware a powerful tool. 

The progression of simulations described here appears 
effective. First, a single narrow slit familiarized stu- 
dents with the simulation and showed more was needed 
for interference. Then, a double slit simulation showed 
Young's experiment in a new light, with visible waves and 
easy measurements. Next, a triple slit experiment made 
students reconsider the constructive and destructive in- 
terference conditions; one held but the other fell. Finally, 
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55% 



25% 



New UWash 



70% 



45% 



My diffraction simulation showed that similar effects could 
a rise from one wide slit. 

While the sample size here was small, the results are 
-encouraging. I hope eventually to try this approach with 



60%. 



Only 15 students used the finial version of the labora- 
tory, so the only conclusion I can draw with confidence 
is that students enjoyed it; I collected anonymous feed- 
back, and it was overwhelmingly positive. Still, for a 
rough gauge of the method's effectiveness, I quizzed the 
students with questions from the University of Washing- 
ton study. Scores for point C are in Table [IJ 

They did well on all parts of the quiz, relative to others. 
60% of the students answered correctly for point C. 80% 
correctly identified the behavior at both points A and 



more students in more classes. 

It would be difficult to experimentally mimic some of 
the simulations with instructional laboratory equipment. 
The simulation could also easily demonstrate more ad- 
vanced topics with minimal modification. The permittiv- 
ity, permeability, and conductivity are definable at every 
point in the simulation domain. Changing those param- 
eters would allow the simulation to demonstrate waveg- 
uides, reflection due to changes in dielectric constants, 
etc. Simulations of a modified two slit experiment — 
where the upper and lower halves of the domain have 
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different indices of refraction — could illustrate more con- 
ceptual topics, like the difference between physical and 
optical path length.^ 
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